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Empirical sizing guidelines such as tail volume coefficients have long been used in the 
early aircraft design phases for sizing stabilizers, resulting in conservatively stable aircraft. 
While successful, this results in increased empty weight, reduced performance, and greater 
procurement and operational cost relative to an aircraft with optimally sized surfaces. 
Including flight dynamics in the conceptual design process allows the design to move away 
from empirical methods while implementing modern control techniques. A challenge of 
flight dynamics and control is the numerous design variables, which are changing fluidly 
throughout the conceptual design process, required to evaluate the system response to some 
disturbance. This research focuses on addressing that challenge not by implementing higher 
order tools, such as computational fluid dynamics, but instead by linking the lower order 
tools typically used within the conceptual design process so each discipline feeds into the 
other. In this research, flight dynamics and control was incorporated into the conceptual 
design process along with the traditional disciplines of vehicle sizing, weight estimation, 
aerodynamics, and performance. For the controller, a linear quadratic regulator structure 
with constant gains has been specified to reduce the user input. Coupling all the disciplines 
in the conceptual design phase allows the aircraft designer to explore larger design spaces 
where stabilizers are sized according to dynamic response constraints rather than historical 
static margin and volume coefficient guidelines. 


Nomenclature 


C D 

= drag coefficient 

deg 

= degrees 

e 

= state vector error 

HT 

= horizontal tail 

i 

= imaginary number 

LQR 

= linear quadratic regulator 

MAC 

= mean aerodynamic chord 

PI 

= Performance Index 

P, q, r 

= x, y, z-body component rotation rates 

RMS 

= root-mean-square 

S 

= wing reference area 

Sh 

= horizontal tail reference area 

S 

= Laplace variable 

s 

= seconds 
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T 

= matrix transpose 

U 

= steady-state velocity in body-fixed stability axis 

u 

= control input vector 

U , V, w 

= x, y, z-body components of velocity 

VT 

= vertical tail 

X 

= state vector 

z 

= performance output vector 

S 

= control surface deflection angle 

a 

= root-mean-square deviation in continuous turbulence 

T 

= time constant 

(j), 0 , ^ 

= roll, pitch, and yaw Euler angles 


I. Introduction 

A S stated by Roskam, “Rather significant improvements in airplane payload-range performance are possible 
by reducing tail areas and relaxing the requirement for inherent static and dynamic stability.” 1 It is 
well known that decreasing the static stability of an airplane will result in overall performance increases. 
Relaxed static stability (RSS), the intentional reduction of static margin, has been used in modern military 
aircraft along with commercial transport aircraft such as the McDonald Douglas MD-11. 2 However, even 
though the benefits are well known, stability and control is rarely applied in the conceptual design process. 
As the demand to reduce fuel burn, emissions, and noise increases to meet the NASA Fixed Wing Project 
next generation aircraft goals, it becomes imperative to incorporate more disciplines into the conceptual 
design process beyond the traditional sizing, weight estimation, aerodynamics, and performance. 

For subsonic flight, the total aircraft parasite drag consists of mostly skin-friction and pressure drag. 3 
With fairings and surface blending, much of the pressure drag can be reduced leaving skin-friction drag as 
the largest contributor to parasite drag. Effectively, if the wetted area of the aircraft can be minimized, the 
skin-friction drag will be minimized resulting in decreased parasite drag. Numerous designs have reduced this 
wetted area by using tailless configurations such as the Convair F-106, Convair B-58, Messerschmitt 163B, 
and the Northrop Grumman B-2. 4,5 All of these configurations have low aspect ratios and thus increased 
induced drag. It would be beneficial if both the induced and parasite drag could be reduced simultaneously. 

A statically stable, conventional configuration requires a download on the horizontal tail for trim, increas- 
ing the required main wing lift coefficient resulting in greater induced drag. On the other hand, positive 
lift produced by the horizontal tail also results in greater than minimal induced drag due to the horizontal 
stabilizer span being less than the main wing. Any lift from the stabilizer results in a total induced drag 
penalty and a reduction in the total aerodynamic efficiency of the aircraft. For a conventional configuration, 
eliminating lift produced by the empennage minimizes aircraft total induced drag. Reducing the aircraft 
static margin is a way to reduce the induced drag on the empennage but requires some form of stability 
augmentation. 6-8 According to Raymer, “[A] modern and sophisticated aft-tail aircraft is designed to a 
slight level of instability so that it normally flies with an upload, not a download on its tail. This is the 
very reason that computerized flight control systems with artificial stability were developed and put into 
production.” 3 

Classical aircraft design uses volume coefficients when sizing an aircraft’s empennage, which tends to 
produce conservative estimates for the stabilizer areas. 9 As a result, the surface area for the horizontal and 
vertical stabilizers exceed the necessary area for adequate dynamic stability. Additionally, this fails to take 
into account the benefits of augmented stability provided by the active control systems in modern aircraft. 
As mentioned previously, a large majority of parasite drag results from skin friction over the wetted area 
of the aircraft. Augmenting the aircraft’s stability through active control will allow the stabilizers to be 
reduced in size thus reducing the total wetted area of the aircraft. In doing so, not only will the induced 
drag be reduced from the reduced tail load but the total parasite drag will be reduced from the decrease in 
wetted area. 

Conceptual design typically focuses on the interaction between the disciplines of aerodynamics, sizing, 
weight estimation, and performance. 10 Any inclusion of control surface design during conceptual design is 
often limited to estimating control surface sizes from historical data while assuming control effectiveness 
is proportional the force and the moment arm. 9,11 Often times, flight dynamics and control (FDC) and 
handling qualities are examined after the aircraft geometry and structural properties have been defined, 
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which inevitably leads to sub-optimal designs or even configurations that have deficient flying qualities. 
FDC is incredibly important when it comes to the overall safety and certification of an aircraft, 9, 12 and 
deficient handling qualities lead to reduced aircraft performance, large cost increases, and delays as the 
configuration must be re-evaluated or redesigned. 

Collaboration between traditional conceptual design disciplines and flight dynamics and controls is es- 
sential when designing for relaxed static stability (RSS). By designing for RSS, aircraft performance is 
increased through reduction of wetted area drag, trim drag, and total weight. 9, n ’ 13 For a transport aircraft 
of conventional stability margins, the horizontal tail accounts for 20-30% of the aircraft-lifting surface and 
approximately 2% of the aircraft empty weight. 11 Any reduction in size of the horizontal stabilizer from 
RSS will provide a significant benefit in reduced drag and aircraft gross weight. However, relaxation of 
stability margin gives a detrimental effect to the aircraft’s handling qualities that must be compensated or 
augmented. 11 It is this correlation between performance gains using RSS and degraded handling qualities 
that make it essential to incorporate flight dynamics and control into the conceptual design phase, especially 
when using a multidisciplinary design optimizer. 

Perez, Liu, and Behdinan in Refs. 9-11 incorporate FDC into a multidisciplinary optimization (MDO) 
including handling qualities, as specified in MIL-F-8785C, on a Boeing 737 class aircraft. The optimizer 
recognized the benefit of RSS by reducing the static margin, moving the wing apex location, shifting the 
CG location, and reducing the horizontal stabilizer (down 28%) and control surface areas. By integrating 
active control into the optimization loop, adequate handling qualities were ensured with feedback control 
to augment the stability. Additionally, reduced control deflections were necessary for trim decreasing trim 
drag. Perez et al. used the decoupled equations of motion with a short period approximation and single-input 
single-output (SISO) to control the dutch roll and short period modes by selecting appropriate gains. 

Morris et al. uses the method of linear matrix inequalities (LMI) to place constraints on the maximum 
actuator deflection, actuator rate, and pole placement limitations. 14 This method relies heavily upon the 
work of Boyd 15 and Kaminer 16 to place constraints on the static feedback gain matrix, K , to obtain desired 
handling qualities. Morris expands his work in Ref. 17 by translating the MIL-STD-1797A guidelines into 
state variance constraints to be used in the development of a state feedback control law using optimal control. 

References 18-21 describe the SirnSAC project using the CEASIOM software which takes a different 
approach than the methods described by Perez et al. and Morris et al. The SirnSAC project uses higher 
order tools, such as computational fluid dynamics (CFD), to iterate a conceptual design. The project has 
had good success showing the benefits of relaxed static stability, but the higher order tools reduce the ability 
to explore a large design space with numerous varying geometric parameters. The optimization time using 
CEASIOM is on the order of weeks instead of the much faster methods described by Perez and Morris and 
the methods presented in this work. 

The focus of this research is to provide a performance gain by reducing the total aircraft drag and weight 
by allowing relaxed static stability and incorporating flight dynamics and control into the conceptual design 
process. By reducing aircraft drag, the overall range and endurance of the aircraft will increase resulting in 
greater operational capability and reduced fuel burn. 

This paper discusses the development of several key components in the MDO framework including flight 
dynamics and control, stability augmentation, atmospheric disturbances, and coupling with a vortex lattice 
aerodynamic code. Fully coupled equations of motion were developed and implemented into a stability 
augmentation system (SAS) that utilizes a linear quadratic regulator where a focus has been placed on 
reducing the user input. Athena Vortex Lattice (AVL) 22 was used to calculate the stability derivatives of a 
Cessna 182 like model to be used in the coupled equations of motion. Unlike much of the work previously 
discussed, state variable transient response to perturbations and atmospheric disturbances were used as a 
measure of handling qualities instead of the more traditional MIL-F-8785C or MIL-STD-1797B military 
specifications. As the configurations and control systems become more complex, the traditional decoupled 
longitudinal and lateral modes may no longer be present, which poses a problem when using pole placement 
techniques. This research eliminates the necessity of identifying the dynamics modes and measuring the 
handling qualities in a qualitative manner and quantitatively evaluates the transient response according to 
the SAE-AS94900 23 specification. Results from varying the horizontal and vertical tail volume coefficients 
are presented. 
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II. Methodology 


A. System Dynamics 

The fully coupled equations of motion were derived and linearized about a steady state condition resulting 
in a state-space representation of the perturbation equations. State feedback was chosen for the controller 
structure due to the guaranteed closed-loop stability properties of the controller which work well in a diverse 
design space. Although it is rare to have all the states available, the focus of this research was not to design 
a robust controller but rather to incorporate active control into the conceptual design space that would give 
a stable solution. 

The perturbation equations were derived assuming a symmetric geometry with zero steady sideslip. 
The steady-state thrust terms were solved for explicitly and substituted into the state-space form of the 
perturbation equations. The implicit, state-space form of the fully coupled perturbation equations is given 

by 

l\l 
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Ex = Ax + 

b : Bg 






where B g is the gust input matrix. The state and control vectors are 
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with the full definition of the matrices in Eq. 1 given in the Appendix. Inverting the generalized inertial 
matrix, E, one obtains the standard, explicit, state-space model 


x = E~ l A + E 


-l 


b : B n 


= Ax 


b : B n 


( 3 ) 


A linear actuator model was added to include the actuator dynamics in the system dynamics. The 
actuator transfer function for each control surface was modeled by a simple-lag filter given by 


M. 

u(s) 


1 

TS + 1 


( 4 ) 


where r is the time constant of the filter. The time constant was chosen to be the same for all three actuators 
and was selected as r = 1/(20. 2) s. 6 The perturbation equations described by Eqs. 1-3 were augmented to 
include the actuator dynamics while neglecting all unsteady terms, thrust terms (except Ct X u ), steady-state 
roll ( R ), and negligibly small terms. AVL was used for the aerodynamic model which does not include thrust 
in the analysis and therefore those terms had to be neglected. The velocity dependent thrust term, Ct Xu , 
was included in the perturbation equations as it would be the largest in magnitude for the configuration 
studied in this research. It was approximated using Ref. 24 to be 


C Txu =~ 3 (C D ~C L sind) 


( 5 ) 


The augmented perturbation equations are presented in compact form in Eq. 6 with the matrix definitions 
given in the Appendix. Again, the hat used in the Appendix equations indicates that a matrix has not be 
premultiplied by the inverted generalized inertial matrix, E aug . 



The actuator input vector, u, is defined as 


B u 



u = 


U e 



( 6 ) 

( 7 ) 
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B. Optimal Control by Linear Quadratic Regulator 

The complete dynamic system of an aircraft incorporating active control is complex requiring numerous 
loop closures to provide adequate closed-loop system response. Classical control relies on the iterative 
selection of gains to achieve the closed-loop system stability, but there is no guarantee the gains chosen will 
be optimum. Modern control theory takes advantage of current computing power where numerous linear 
equations can be solved simultaneously to obtain a set of gains that minimizes a chosen performance index 
(PI). It is in the selection of the performance index that the true engineering of the controller occurs. 

A linear quadratic regulator (LQR) can be used to simultaneously close all the loops in a linear, time- 
invariant (LTI), multi-input multi-output (MIMO) system. With the closure of all the loops, the gains are 
solved simultaneously negating the need for successive loop closure as required in classical control theory. 
Extremely versatile, the LQR is capable of using performance indices with state and control weighting, 
time weighting, and derivative weighting of the state in both state feedback and output feedback control 
structures. Without any restrictions on the gains, other than closed-loop stability is required, the LQR may 
choose to zero a gain, thus leaving a loop unclosed. Additionally, a compensator may be used in the form of 
filters, integral and derivative controllers. This flexibility makes it an excellent tool for finding optimal gains 
for a controller in an LTI system. 

As a regulator, any non-zero states are driven to zero in such a way that a chosen performance index is 
minimized. By driving all the states to zero the system is returned to the steady-state condition with the 
least amount of cost, thus being an optimal solution. This is ideal for a stability augmentation system where 
any deviation from the steady state is undesired. 

A set of nonlinear equations of motion can be linearized about a steady-state condition and represented 
in state-space form as given by 

x = Ax + Bu (8) 

where both x and u are functions of time. The state vector is a vector of perturbations from the steady state 
condition that the regulator drives to zero. With the state-feedback control law, 

u = —Kx (9) 

the closed loop system takes the form 

x = {A — BK) x = A c x (10) 

A performance output can be defined as 5 6 

z = Hx (11) 

where 2 is a combination of states. For a regulator, z can be set to the error of a specific state such as a 
non-zero pitch rate that should be driven to zero. 

The LQR finds the optimal gains through the minimization of a performance index that integrates the 
values of both the state and control vectors over time. The impact of each state and control is weighted in 
the performance index through the use of weighting matrices. The standard performance index for the LQR 
is 6 

1 f°° 

J = - / (x T Qx + u t Ru) dt (12) 

2 J o 

with Q > 0, R > 0. The performance output, z, can be incorporated into the performance index such that 
Q = H t H. 

A limitation of the linear quadratic method is the n x n entities that must be chosen in the weighting 
matrix Q where the values may not correspond to a performance objective due to an observability require- 
ment, initially presented by Kalman 25 and discussed by Stevens and Lewis. 6 This results in a trial-and-error 
method of selection for Q where the entries are varied until an acceptable transient response is obtained. 
This method of design is highly undesirable. 

Eliminating the restriction of the observability requirement in the selection of Q allows for the entries to 
be chosen strictly on desired performance objectives. With specified performance objectives, the structure of 
the PI and the number of entities to be chosen for the weighting matrices can be reduced and the closed-loop 
response is dependent on the design of the performance index. 

A strength of the linear quadratic method is the flexibility in the selection of the performance index 
structure. The standard PI given in Eq. 12 can be modified by adding time and derivative weighting of the 
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states. A benefit of the time-weighted performance index is that it can satisfy the observability requirement 
freeing the selection of the weighting matrices. The standard PI only lightly penalizes small errors due to 
a slow pole(s) with small residue resulting in the time to reach a steady state condition to be rather large. 
The time-weighted performance index heavily penalizes errors that occur late in the response and as a result 
suppresses the effect of a slow pole and lightly damped settling behavior. 6 

The derivative of the state can sometimes be a more accurate representation of the workload on a control 
system and should be weighted in the performance index rather than the state itself. For example, the rate 
of change of the elevator is a more accurate representation of required actuator power than the deflection 
angle itself. 

The standard PI for the linear quadratic regulator was modified to include the state error and derivative 
weighting of the states given by 




e + u t Ru + z T Wz ) dt 


(13) 


where R has the same definition as in the standard performance index, Q has been set equal to zero, W is 
a new weighting matrix on the states’ rate of change. The state error, e, is defined as 


e = z = Hx 


where P = H T H. 

The solution to Eq. 13 can be shown to be 

J = 1 tr (PkX) 

subject to the nested Lyapunov equations given by 

0 = go — AjPo + PqA c + P 
0 = gi = AjPi + P\A C + P 0 

0 = Qk - 1 = Aj Pk-1 + Pk-\A C + Pfc — 2 
0 = g k = A T c P k + P k A c + fc!P fc -i + K t RK + A T C H T WHA C 


(14) 


(15) 


(16) 


with X defined as 

X = x(0)x T (0) = H~ 1 z(0)z t (0) (P- 1 ) T (17) 

The solution to this minimization problem is dependent upon the initial conditions, x(0). A simple 
solution for eliminating the dependence on the initial state is to average the performance of a set of linearly 
independent initial conditions. This is equivalent to a random variable uniformly distributed on the surface 
of a unit sphere of dimension equal to the length of the state vector. 26 In essence, instead of minimizing the 
performance index in Eq. 15, the expected value of J, E{ J}, is minimized such that 

E{J} = ±E{x t (0)Px(0)} = \tr{PX) (18) 


where the symmetric n x n matrix 


X = P{x(0)x T (0)} 


(19) 


is the initial autocorrelation of the state. For the regulator problem, it is practical to set X = I since it is 
desired to drive arbitrary nonzero states to zero. 6 

A emphasis was placed on reducing the number of entities that must be chosen in the weighting matrices. 
As a regulator with state feedback, any nonzero state was considered an error. The performance output, z, 
was selected as 

r n T 


uj3apqr(j)6'il)5 e 


(20) 


with the performance output matrix H defined as 


H = 


180 


diag | 


— A A 1 1 1 1 

180 u u 


11111 


(21) 
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The performance output was chosen to give one degree of error equal weighting as one foot per second in the 
performance index. This was chosen to maintain the generality of the controller while reducing the number 
of entities to be chosen in the performance index. 

The derivative-weighting matrix W was chosen as 

W = diag{o 0000000011 lj (22) 


to place a small penalty on high control surface deflection rates. 
The controls- weighting matrix R was selected as 


R = 0.1 


1 0 
0 1 
0 0 


0 

0 

1 


(23) 


where smaller values of R gives greater authority to the control inputs. However, values less than those given 
in Eq. 23 produce transient responses negligibly different than results using Eq. 23. 


C. Atmospheric Disturbances 


To stress the control system, atmospheric disturbances 
were modeled and used to test the response of the system to 
deviations from the steady-state condition. Including these 
disturbances adds validity to the model and eliminates un- 
realistic configurations. Two atmospheric disturbance mod- 
els were used as suggested in military standard MIL-STD- 
1797A: 27 a continuous turbulence model in the frequency 
domain and a discrete gust model in the time domain. 


1. Continuous Turbulence 

The von Karman continuous turbulence model was used 
as specified in MIL-STD-1797A. The spectra are of the form 


</S,( n ) = T, 


<t>va(ty = 


2 


i + (1.339 L u ny 


n 5/6 


1 + f (1.339i„H) 2 


1 + (1.339L„H) : 


11/6 


(24) 


Table 1. RMS gust intensities, fps. 


Altitude (ft) 

Light 

Moderate 

Severe 

500 

6.6 

8.6 

15.6 

1,750 

6.9 

9.6 

17.6 

3,750 

7.4 

10.6 

23.0 

7,500 

6.7 

10.1 

23.6 

15,000 

4.6 

8.0 

22.1 

25,000 

2.7 

6.6 

20.0 

35,000 

0.4 

5.0 

16.0 

45,000 

0 

4.2 

15.1 

55,000 

0 

2.7 

12.1 

65,000 

0 

0 

7.9 

75,000 

0 

0 

6.2 

85,000 

0 

0 

5.1 


KniP) = Cl 


i + | (1.339 L w ny 


i + (1.339 L w ny 


1 11/6 


where L n is the length scale and a n is the root-mean-square intensity of the continuous turbulence. The 
spatial frequency, f 2, used in the von Karman form of the spectra is related to the temporal frequency, oj, by 
f l = co/U. The spectra are converted to functions of u> using the simple relationship 




n = u.v.w 


(25) 


The mean square of the output can be found be integrating the square of the system transfer function iui 
times the power spectral density of the continuous response for all frequencies as shown in Eq. 26. 


= 2 


\G(iu)\ 2 (j> n (ui)dui, 


n = u.v.w 


(26) 


For isotropic atmospheric turbulence, a requirement of both the von Karman and Dryden forms, the turbu- 
lence length scales and mean-square intensities are the same for the three velocity components and is shown 
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by 


2 2 2 
07 = at = cC 


Lit. — Li 1 — Li j 


(27) 


For clean air turbulence above 2,500 feet, the length scale of the turbulence is 2,500 feet in the von Karman 
form and 1,750 feet for the Dryden form which is used in the discrete gust model described in the next 
section. 

The RMS intensities to be used in Eq. 24 were specified in Table 4 in Section 3. 2. 1.5.1 of SAE-AS94900. 23 
The table gives the RMS intensities for numerous altitudes with varying probability of exceedance. The 
intensities corresponding to light, moderate, and severe with probabilities of exceedance of 10~ 2 , 10 -3 , and 
10“ 5 are given in Table 1. Altitudes used between data points were linearly interpolated to allow for any 
desired altitude. 


2. Discrete Gust 

The discrete gust model has the “1-cosine” profile defined by Eq. 28 and illustrated in Fig. 1 


v = 0 , 


x < 0 


V m ( ttx . 

v = — II — cos — ) , 0 < X < d m 


(28) 


v = 0, 


X > dr, 


where V m is the magnitude of the gust and d m is one-half the total gust length. The gust length was chosen 
to provide the maximum system excitation by tuning the gust to the system natural frequencies. 23 The 
magnitude of the gust was chosen according to MIL-STD-1797A using Fig. 2. 

The Dryden turbulence scales, and the root- 
mean-square intensities from Table 1, are used 
in Fig. 2. SAE-A94900 specifies the maximum 
discrete gust to be a single full wave with a peak 
amplitude of 60 fps. Figure 2 is for light to 
moderate turbulence where the probability of 
exceedance is 10~ 2 . The user of the tool has the 
option to use either the calculated gust magni- 
tude taken from Fig. 2, or use the severe turbu- 
lence case of 60 fps. 




NORMALIZED DISCRETE GUST LENGTH, ^ \ d z 

Lu Ly Lw 

Figure 2. Normalized discrete gust for determining gust magnitude. 27 ’ 29 
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D. Dynamic Performance Constraints 

Reducing horizontal and vertical stabilizer areas results in degraded flying qualities of any configuration. 
Damping of the dynamic modes is directly affected by those stabilizer areas. The ability to trim any moments 
also reduces with the empennage size due to reduced control effector area. Flying quality enhancements must 
be provided through active control systems that reduce or eliminate undesired residual oscillations. The 
control system workload increases with decreasing empennage area resulting in a point where the control 
system cannot provide adequate flying qualities. 

A 20 degree maximum control deflection angle was chosen to limit the authority of the control system. 
Twenty degrees was chosen to allow a margin of safety in addition to avoiding the nonlinear effects of large 
control surface deflections. The total deflection angle was calculated as the sum of the steady-state deflection 
and any perturbation deflections. This was checked during the discrete gust simulation discussed previously. 

Military specification SAE-AS94 9 00 23 specifies performance requirements on the transient response of 
the attitude angles that allows for quantitative evaluation of the flying qualities without having to identify 
specific dynamic modes; this is advantageous because the traditional dynamics modes may not even exist 
due to the complex dynamics provided by active control. Assessing the transient response allows for the 
flying qualities to be evaluated for any geometric configuration with any active control system. 

As specified in SAE-AS94900, the root- mean-square deviations in pitch attitude must be less than or 
equal to five degrees in a continuous turbulence field. The control system must also be capable of returning 
the pitch attitude to within plus or minus 0.5 degrees of the steady-state condition, resulting from a five 
degree perturbation in pitch, within five seconds for aircraft in classes I— III, defined in MIL-STD-1797A. 
The roll attitude requirements are similar to the pitch attitude requirements with the exception of a ten 
degree root-mean-square deviation in roll requirement in continuous turbulence, and a static accuracy of one 
degrees within five seconds from a roll perturbation. In continuous turbulence, the heading angle, i/>, must 
have a root-mean-square heading deviation of less than or equal to five degrees. 23 

While climbing up to a maximum rate of 2000 feet per minute, the control system must be capable 
of leveling off and achieving a static airspeed accuracy of plus or minus 10 knots or 2% of the reference 
airspeed, whichever is greater. This accuracy must be achieved within 30 seconds of engaging the airspeed 
hold. Any residual oscillations within the static accuracy margin must have a period of oscillation greater 
than 20 seconds. This requirement was modeled as a small perturbation in the pitch angle, 6. 

Table 2. Summary of dynamic performance constraints. 


Description 

Test Condition 

Constraint 

20 deg max control deflection 

pert, and gust 

never exceed 

5 deg pitch 

perturbation 

± 0.5 deg in < 5 s 

5 deg roll 

perturbation 

± 1.0 deg in < 5 s 

airspeed hold 

perturbation 

± 10 kts or 2% < 30 s 

pitch deviation 

cont. turbulence 

erg < 5 deg 

roll deviation 

cont. turbulence 

(70 < 10 deg 

heading deviation 

cont. turbulence 

(70, < 5 deg 


E. Flight Condition 

Only one flight condition was modeled for the development of this tool and was used in the initial results 
presented in Section IV. The cruise flight condition was implemented in the development of the tool due to 
the availability of aerodynamic data given by Napolitano. 30 This aerodynamic data was used to verify the 
derived equations of motion by matching the roots of the state matrix to the roots of the transfer function 
characteristic equation derived in Ref. 30. Additionally, data on the flight condition such as altitude, Mach 
number, true airspeed, steady-state angle of attack, and dynamic pressure were given. Finally, the mass 
properties were given in Napolitano, a key component to accurately evaluating the dynamics of any aircraft. 
Data taken from Napolitano are summarized in Table 3. 

The cruise condition was chosen only for the development of the tool. It is well acknowledged that this 
flight condition is not a key flight condition for sizing of the stabilizers and other flight conditions must be 
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added to increase the validity of the analysis. 


Table 3. Cessna 182 cruise flight condition. 30 


Altitude 

(ft) 

u 

(ft/s) 

M Dynamic Pressure 

(-) (IWft 2 ) 

a 

(deg) 

m 

(ib) 

lx 

(si - ft 2 ) 

Iy 

(si - ft 2 ) 

Iz 

(si - ft 2 ) 

5,000 

220.1 

0.201 49.6 

0 

2,650 

948 

1,346 

1,967 


III. Analysis 


A. Aerodynamic Modeling 


Athena Vortex Lattice (AVL), a code 
maintained by Mark Drela at the Massachusetts 
Institute of Technology (MIT), was used to eval- 
uate the aerodynamics of each configuration. 

This tool was chosen due to the ease of availabil- 
ity, previous user experience, ability to calculate 
stability derivatives while trimmed in both lat- 
eral and longitudinal directions simultaneously, 
and the low computational cost. The stability 
derivatives were taken directly from the AVL 
output and input into the equations of motion 
described in Section II. 

The Cessna 182T Skylane geometry was 
modeled using 3-view sketches from Refs. 31 
and 32. Only the lifting surfaces were modeled 
in AVL while the non-lifting surfaces were ne- 
glected to maintain simplicity of the geometry. The control surfaces were accurately modeled and then 
modified so the control surfaces were defined as a constant percent chord. This reduces the number of 
geometry parameters in the tool while the effects on the analysis were insignificant. 

Figure 3 shows the Cessna Skylane model pictured in AVL’s graphical user interface (GUI) with the 
non-lifting bodies removed. The center of gravity (CG) position was placed 8.851 feet measured from the 
propeller spinner as shown in the Pilot’s Operating Handbook 31 and given by Napolitano. 30 



Figure 3. AVL model of Cessna 182 Skylane. 


B. Modification of the Geometry 


In this initial study, the only geometrical parameters varied were the horizontal and vertical tail volume 
coefficients. The vertical tail strake, shown in Fig. 3, was neglected to maintain the trapezoidal geometry 
parameters. To reduce the number of variables, the horizontal and vertical tail moment arms, aspect ratios, 
quarter-chord sweeps, dihedral angles, taper ratios, and airfoils were unchanged from the baseline model. 
Regardless of the change in volume coefficient, the root quarter-chord coordinate locations were identical. 
Effectively, changing the volume coefficient with these restrictions results only in a change of reference area 
and span. 

The horizontal tail volume coefficient is defined as 


_ ShIh 
Sc 


(29) 


Traditionally, Ih is defined as the length between the wing aerodynamic center and the stabilizer aerodynamic 
center. For simplicity due to the current fixed CG position in this analysis, the length 1^ used in Eq. 29 
was measured from the CG position to the root quarter-chord point on the stabilizer. Of note, the volume 
coefficient for the vertical stabilizer is the same equation as Eq. 29 with the subscript H replaced by subscript 
V and mean aerodynamic chord, c, replaced by the span, b. 

The tail volume coefficients were varied from the Cessna 182 model values of Ch — 0.7 and Cy = 0.28 
to 25% of the baseline over ten equally spaced steps. This resulted in 100 geometries evaluated by the tool, 
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each with differing stabilizer areas and spans, taking nearly 50 minutes of computation time. As the tool 
was still in development, computation time was sacrificed to provide extra data to the user, ensuring it was 
performing as expected. 



Figure 4. Cessna 182 geometry with modified horizontal 
tail as percentage of baseline volume coefficient. 


Figure 4 shows a drawing from the 
top view of the AVL geometry with 
varying horizontal tail volume coeffi- 
cients. The coefficients shown are 0.70, 
0.52, 0.35, and 0.17 which correspond to 
100%, 75%, 50%, and 25% of the base- 
line. The centerline of the fuselage is in- 
dicated by the solid line drawn through 
the symmetric plane. 

Any experienced aircraft designer 
would notice the small horizontal tail 
area for the 25% baseline volume coef- 
ficient and doubt its ability to achieve 
adequate handling qualities, much less 
trim at takeoff conditions. As discussed 
later in Section IV, all cases passed the 
dynamic performance checks for ade- 
quate handling qualities. However, fo- 
cusing only on the cruise flight condi- 
tion is inadequate for accurately sizing 
the stabilizers, and more constraining 
flight conditions must be added to the 
analysis. 


C. Drag Buildup 


To capture all of the drag benefit of reducing the stabilizer area, and including active control in the design 
process, the reduction in parasitic drag must be calculated from the reduced wetted area. A zero-lift drag 
coefficient in the cruise condition was given by Napolitano as Cd 0 = 0.027. This is the total parasitic drag for 
the Cessna 182 Skylane which was modeled using AVL. Reducing the horizontal and vertical stabilizer areas 
will reduce the total parasitic skin-friction drag. To model this, the contribution of the baseline empennage 
to the total parasitic drag was estimated using Raymer’s component buildup method where each component 
is denoted by the subscript c. 3 The subsonic parasitic drag coefficient can be approximated by 


(C Do ) 


subsonic 


E Cf c FF c Q c S wetc 
s 


(30) 


where C/ is the skin- friction drag coefficient, FF is a form factor that estimates the pressure drag due to 
viscous separation, Q is an interference factor, and S we t is the estimated wetted area of the component. The 
miscellaneous, leakage, and protuberance drag terms were ignored. The empennage contribution to Cd 0 was 
subtracted from the baseline configuration Cn 0 of 0.027 leaving only the wing-body parasitic drag. This 
allows the stabilizer areas to be varied and the effect on the parasitic drag to be captured. Hence, reduced 
stabilizer areas result in reduced parasitic drag. 


IV. Results and Discussion 

Figure 5 shows the total drag coefficient as a function of the horizontal and vertical tail volume coefficients. 
As expected, the total drag decreased as the tail volume coefficients were reduced; the total drag was sensitive 
to both volume coefficients, but had a greater sensitivity to the vertical tail volume coefficient. Figure 6(a) 
shows the reduction in drag that resulted from decreasing the static margin. The drag benefit seen in 
Figs. 5 and 6(a) comes from the reduction in wetted area and the resulting decrease in parasitic drag. When 
trimmed, induced drag for this geometry actually increased with the reduction of static margin as shown in 
Fig. 6(b). This explains why the total drag was more sensitive to varying the vertical tail volume coefficient 
as mentioned previously. The resulting increase in induced drag comes not from the reduction of static 
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margin but the reduction in the span, and increase in lift coefficient, of the horizontal stabilizer. 
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Figure 5. Drag coefficient for varying tail volume coefficients. 

The elevator deflection angle, shown in Fig. 7, increased in magnitude as the static margin was decreased. 
The Cessna 182 baseline geometry required zero elevator deflection to trim at the cruise condition which 
means the downwash from the main wing provided sufficient induced angle on the horizontal stabilizer to 
trim. Because the location of the wing and center of gravity were fixed, the horizontal tail has to generate the 
same amount of moment to trim as the baseline case. As the horizontal tail area and span were decreased, 
the down lift on the stabilizer decreased requiring the elevator to deflect, increasing the magnitude of the 
horizontal tail lift coefficient, to achieve trim. As a result, a net increase in configuration induced drag was 
realized. This indicates that simply reducing the static margin does not decrease induced drag but that a 
system level design, with static margin being allowed to decrease, must be performed in order to achieve the 
induced drag benefits of relaxed static stability. 



Figure 6. Sensitivity of total and induced drag coefficients to static margin. 


The accuracy of the aerodynamic model should also be considered when designing for relaxed static 
stability. AVL predicted a baseline static margin of nearly 30% which was much greater than the actual 
static margin. Both the engine and fuselage, neglected in the AVL model, will shift the neutral point forward 
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decreasing the static margin. This modeling error causes the geometry to be overly stable in the analysis of 
the dynamic response. 

Initially, using moderate atmospheric disturbances with a probability of exeedance level of 10 — 3 , none of 
the configurations passed the heading hold turbulence check specified in SAE-AS94900. This resulted from 
the active control system not sufficiently driving the rudder in continuous turbulence. A great feature of 
the LQR is the flexibility in allowing for adjustments in the weighting matrices if a desired performance is 
not achieved. Due to the consistent failures in the heading hold check, the ip weighting term in Eq. 14 was 
increased by a factor of 10, placing heavier penalties on non-zero values of yaw. This simple adjustment 
enabled all configurations to achieve adequate dynamic performance in the cruise flight condition. It should 
be noted, again, that the cruise flight condition alone does not adequately stress the control system for sizing 
the stabilizer surfaces. Additional flight conditions must be added to ensure proper sizing of the stabilizer 
surfaces. 

Figures 8-12 show the baseline con- 
figuration’s dynamic response with the 
heavier yaw weighting to perturbations, 
continuous gust fields, and discrete 
gusts. These disturbances were used to 
check the dynamic performance of the 
active control system as summarized in 
Table 2. The airspeed hold and lon- 
gitudinal gust responses were omitted 
for brevity due to their similarity to the 
pitch hold and vertical gust responses. 

Time histories of the transient re- 
sponse to the perturbations specified 
by SAE-AS94900 are shown in Figs. 8 
and 9. The figures were separated into 
longitudinal and lateral states, show- 
ing both the short time response — less 
than ten seconds — and the long time 
response — total response over 100 sec- 
onds. All control surface states shown 

in Figs. 8, 9, 11, and 12 are the total Fi § ure 7 ' Elevator deflection angle versus static margin. 

deflection angles, the sum of the perturbation and the steady-state deflection angle. Clearly visible in the 
open-loop response of Fig. 8 is the phugoid mode which has clearly been attenuated with little control 
power. Figure 9 shows the unstable spiral mode in the open-loop response that has been stabilized through 
a combination of aileron and rudder control surface deflections. 

Power spectral densities of the mean-squared response are plotted in Fig. 10. With a natural frequency 
of 0.19 rad/s, the phugoid sensitivity can clearly be seen in the vertical continuous turbulence responses as 
a spike in the mean-squared response of the open-loop system. Less obvious is the short period mode, which 
had an open-loop natural frequency of 7.66 rad/ s , that can be seen as a slight hump in the angle of attack, 
pitch rate, and pitch angle. The spike in the response of the lateral states of Fig. 10 correspond to the dutch 
roll natural frequency of 5.14 rad/s. The active control system successfully reduced the magnitude of the 
response of the phugoid and dutch roll modes while the short period response in turbulence was unchanged. 

Figures 11 and 12 show the response to a lateral and vertical discrete gust. As with the pitch pertur- 
bations, the phugoid is quite visible in the open-loop response. Due to the model being extremely stable 
at cruise, even in a discrete gust, the required elevator deflection of the closed- loop system was extremely 
small while providing a heavily damped response. The closed-loop lateral gust response, Fig. 11, provides 
increased damping when compared to the open-loop response, as would be expected. Heavier weighting in 
the performance index on the heading angle, ip, required to meet the required heading hold performance 
results in a sacrifice in the roll angle, obtaining improved heading angle response. The heading perturbation 
state of Fig. 11 deviates from steady-state in response to the gust but immediately returns to the steady-state 
value, whereas the roll perturbation state has two full oscillation cycles before returning to the steady state. 



Static Margin, %MAC 
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Figure 8. Pitch hold perturbation check. 
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Continuous Turbulence Response 



Figure 10. Continuous turbulence response. 
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Figure 11. Lateral discrete gust response. 
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Figure 12. Vertical discrete gust response. 
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V. Conclusions 

A multidisciplinary optimization, conceptual design tool methodology was presented along with initial 
results. An LQR controller was developed along with a methodology for selecting the weighting matrices 
that reduces number of elements that must be chosen. Using control system performance requirements from 
SAE-AS94900, dynamic constraints were implemented and tested using perturbation simulations, continuous 
atmospheric turbulence modeling, and discrete gust simulations. The geometry was varied and active control 
used to stabilize the system to allow for a reduction in wetted area and static stability. Initial results showed 
that a heavier weighting on the heading state was required to meet the continuous turbulence performance 
requirements of SAE-AS94900. After the adjustment in heading state weighting, drag reductions were 
achievable through the reduction of the horizontal and vertical stabilizers while maintaining acceptable 
dynamic response. By implementing active control using a linear quadratic regulator, with a simplified 
structure to reduce weighting matrix decisions, the design process, including active control, has been shown 
to be feasible. As indicated in the results, a decrease in static margin does not guarantee a reduction in 
induced drag; to realize the full benefits of relaxed static stability, the entire configuration needs to be given 
the freedom to vary geometric parameters without the traditional restriction on minimum static margin. 

Continued expansion of the tool is necessary to incorporate all the traditional design disciplines to make 
it fully functional. The results indicate the importance of identifying the constraining flight conditions when 
sizing the stabilizers as the cruise flight condition did not heavily load the active control system. 

VI. Future Work 

Additional flight conditions will be added as the cruise condition used here was not a constraining flight 
condition for sizing control effectors. The atmospheric disturbances model will also be enhanced to include 
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the induced rotation rates resulting from each surface entering the disturbance at slightly different times. 
This effect will be large for aircraft with large separations between the main wing and horizontal stabilizer. 
The discrete gust magnitude was calculated using Ref. 28 and had a probability of exceedance of 10 — 2 , a light 
disturbance level. Greater disturbance levels will be added to the normalized discrete gust calculation to 
account for more severe gust levels. These additions will increase the validity of the atmospheric disturbance 
modeling and place a greater workload on the active control system. 

The work presented herein gave an indication of the intentions of this multidisciplinary design tool. A 
sizing and weights estimation tool needs to be added to allow for the variation of the entire geometry in the 
design space. The Flight Optimization Software (FLOPS) code will be attached to the method discussed 
and will perform the sizing, weight estimates, performance estimates, and parasitic drag buildup. A baseline 
FLOPS model of an advanced, 737-800 class, double-bubble fuselage configuration called the D8 will be used. 
A global optimizer will then be allowed to vary numerous geometric parameters with the goal of including 
flight dynamics and control into the conceptual design process, helping to achieve the Aeronautics Research 
Mission Directorate goal of reduced fuel burn. 


Appendix 

The matrices of the fully coupled, perturbation equations are provided below. These matrices are used 
in Eqs. 1, 3, and 6. 
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The symbols used 

in Eqs. 31-39 are defined as 



b 

wingspan 

s 

= wing reference area 

Qoo — 

dynamic pressure 

m 

= aircraft mass 

a = 

steady-state angle of attack 

g 

= gravity constant 

R 

steady-state roll rate 

c 

= mean aerodynamic chord 

ii 

stability axis mass 

T' 

1 XZ 

= stability axis cross- 


moments of inertia 


product of inertia 
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